Search for three-nucleon force effects 
on the longitudinal response function of 4 He 
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A detailed study of the 4 He longitudinal response function Rt,{u),q) is performed at different 
kinematics, with particular emphasis on the role of three-nucleon forces. The effects shown are the 
results of an ab initio calculation where the full four-body continuum dynamics is considered via the 
Lorentz integral transform method. The contributions of the various multipoles to the longitudinal 
response function are analyzed and integral properties of the response are discussed in addition. The 
Argonne Vf 8 nucleon-nucleon interaction and two different three-nucleon force models (Urbana IX, 
Tucson- Melbourne') are used. At lower momentum transfer (q < 200 MeV/c) three-nucleon forces 
play an important role. One even finds a dependence of Rl on the three-nucleon force model itself 
with differences up to 10%. Thus a Rosenbluth separation of the inclusive electron scattering cross 
section of 4 He at low momentum transfer would be of high value in view of a discrimination between 
different three-nucleon force models. 

PACS numbers: 25.30.Fj, 21.45.-v, 27.10.+h, 31.15.xj 
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I. INTRODUCTION 

An aspect of nuclear dynamics that has attracted a lot 
of interest in the last years is the importance of multi- 
nucleon forces and in particular of the three-nucleon force 
(3NF). The nuclear potential has clearly an effective na- 
ture, therefore it is in principle a many-body opera- 
tor. Yet the debate has concentrated for several decades 
mainly on its two-nucleon part. Such a debate has taken 
place among the advocates of three different approaches 
based on meson theory, pure phenomenology and more 
recently effective field theory. Realistic potentials have 
been obtained within the three different frameworks, re- 
lying on fits to thousands of N-N scattering data. As 
is well known such realistic potentials do not explain the 
triton binding energy and thus 3NF are necessary. Today, 
due to effective field theory approaches a new debate is 
taking place regarding 3NF. However, for the determina- 
tion of a realistic three-body potential or to discriminate 
among different models one needs to find A> 3 observ- 
ables that are 3NF sensitive. An important activity in 
this direction has taken place in the last years, with ac- 
curate calculations of bound-state properties of nuclei of 
increasing mass number A [l], Q ■ 

We follow a complementary approach and direct our 
attention, instead, towards electromagnetic reactions in 
the continuum. In fact many years of electron scatter- 
ing experiments have demonstrated the power of electro- 
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nuclear reactions, and in particular of the inelastic ones, 
in providing important information on nuclear dynam- 
ics. The possibility to vary energy ui and momentum q 
transferred by the electron to the nucleus allows one to 
focus on different dynamical aspects. In fact one might 
find regions where the searched three-nucleon effects are 
sizable. The choice of the 4 He target is particularly ap- 
propriate because of the following considerations, i) The 
ratio between the number of triplets and of pairs goes like 
(A-2)/3, therefore doubles from 3 He to 4 He. ii) Theoret- 
ical results on hadronic scattering observables involving 
four nucleons have already shown that three-body effects 
are rather large [1, Q| ■ hi) Since 4 He has quite a large av- 
erage density and its binding energy per particle is simi- 
lar to that of heavier systems, it can serve as a guideline 
to investigate heavier nuclei, iv) Various inclusive 4 He 
(e, e') experiments have already been performed in the 
past Q, where Rosenbluth separations have been car- 
ried out. Due to the low atomic number it is possible to 
study the longitudinal and the transverse responses sep- 
arately, without the ambiguities created by the Coulomb 
distortions affecting heavier systems, v) The Lorentz In- 
tegral Transform (LIT) method @, allows to extend 
investigations beyond the three- and four-body break up 
thresholds. 

In this work we concentrate on the longitudinal re- 
sponse function R^^uj^q) at constant momentum trans- 
fers q < 500 MeV/c. Since the longitudinal response 
Rl is much less sensitive to meson exchange effects than 
the transverse response Rt the use of a simple one-body 
density operator allows to concentrate on the nuclear dy- 
namics generated by the potential. In fact, for low q two- 
body operators in Rl are only of fourth order in effective 
field theory counting (N 3 LO) [8(, and their contribution 
is negligible up to q w 300 MeV/c, see Sec. IVl 
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Besides presenting new results this work gives a more 
detailed analysis of those published in a previous Let- 
ter [§]. The paper is organized as follows. In Sec. Hi] 
we give the definition of Rl and explain the theoreti- 
cal framework that allows to calculate it. In Sec. IIIII we 
show the results for different kinematics and compare our 
results with existing data. In Sec. IIVI we analyze our re- 
sults as obtained from a multiple decomposition of the 
response function. In Sec [V] we discuss integral proper- 
ties of the longitudinal response and compare them with 
some of the results in the literature. Finally conclusions 
are drawn in Sec. IVII 



II. THEORETICAL FRAMEWORK 

In the one-photon-exchange approximation, the inclu- 
sive cross section for electron scattering off a nucleus is 
given in terms of two response functions, i.e. 
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where <jm denotes the Mott cross section, Q' 2 = — q 2 = 
q 2 — lu 2 the squared four momentum transfer with ui 
and q as energy and three-momentum transfers, respec- 
tively, and 9 the electron scattering angle. The longi- 
tudinal and transverse response functions, RL(uJ,q) and 
i?r(w, q), are determined by the transition matrix ele- 
ments of the Fourier transforms of the charge and the 
transverse current density operators. In this work we 
focus on the longitudinal response which is given by 



R L (u,q) = £ | (V f \p(q)\* ) \ 2 S [ej + ±- _ E{ 



Q—U) 



(2) 

where M is the target mass, |\&o//) an d -^o// denote ini- 
tial and final state wave functions and energies, respec- 
tively. The charge density operator p is defined as 



Kl) = | + T k) ex P t ic t ' rfc l ' 



(3) 



where e is the proton charge and r| the isospin third 
component of nucleon k. The ^-function ensures energy 
conservation. 

As it will be clear in Sec. IIIII it is useful to consider the 
charge density operator as decomposed into isoscalar (S) 
and isovector (V) contributions 

p(q) = f Yl exp t* q ' rfc ] + o H r fc exp t iq ' rfe ] 



As(q) +/V(q)- 



(4) 



Each of them can be further decomposed into Coulomb 
multipoles [l(| 



p x ( (l )=inJ2^' X (q)Y l f(qy 

J ft 



(5) 



where the Coulomb multipole operators C^ ,x (q) are de- 
fined, by 



4tt 



dq'p x (q')Y u J (q'), 



(6) 



with X = S, V and Yjf(q) denoting the spherical har- 
monics. 

From Eq. @ it is evident that in principle one needs 
the knowledge of all possible final states excited by 
the electromagnetic probe, including of course states in 
the continuum. Thus, in a straightforward evaluation 
one would have to calculate both bound and contin- 
uum states. The latter constitute the major obstacle 
for a many-body system, since the full many-body scat- 
tering wave functions are not yet accessible for A > 3. 
In the LIT method @, 0] this difficulty is circumvented 
by considering instead of Rl{oj, q) an integral transform 
Cl (<r, q) with a Lorentzian kernel defined for a complex 
parameter a = uji + i 07 by 



Cl (a, q) = dui 



(w - ctr) 2 + a 2 V a ' 9 ' a - q/ 



(7) 



The parameter ai determines the resolution of the trans- 
form and is kept at a constant finite value (07 ^ 0). The 
basic idea of considering Cl lies in the fact that it can 
be evaluated from the norm of a function , which is 



the unique solution of the inhomogeneous equation 
(H-E -a)\*r tq )=p(q)\* ). 



(8) 



Here H denotes the nuclear Hamiltonian. The existence 
of the integral in Eq. implies that ^ has asymptotic 
boundary conditions similar to a bound state. Thus, one 
can apply bound-state techniques for its solution. Here 
we use the effective interaction hyperspherical harmonics 

(EIHH) method [ll 03]. 

The response function RL(t-u,q — const) is then ob- 
tained by inverting the integral transform ([7]). For the 
inversion of the LIT various methods have been de- 
vised [HI, In particular the issue of the inversion 
of the LIT is discussed extensively in Ref . [l5[ • 

Finally we should mention that the expression of the 
charge density in Eq. ([3]) describes point particles. In or- 
der to compare our results with experimental data, after 
inversion the isoscalar and isovector parts of Rl have to 
be multiplied by the proper nucleon form factors 

1 -(1 + tI)->G s e {Q 2 ) + tIG v e {Q 2 ), (9) 

where G E and G\ are the isoscalar and isovector form 
factors 



G S E = l(G p E + G n E ), 



2 V E 
5^ 
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(11) 



For on-shell particles, these form factors depend on the 
squared four momentum transfer Q 2 alone. In principle, 
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FIG. 1: (Color online) Longitudinal response function for 
q = 200 and 100 MeV/c with the AV18 (dashed), AV18+UIX 
(solid) and MT (dashed-dotted) potentials. Data at q ~ 200 
MeV/c from [r7[. 



this is no longer true for the off-shell situation. How- 
ever, in view of the fact that little is known about the 
off-shell continuation and, furthermore, for the moderate 
energy and momentum transfers considered in this work, 
the neglect of such effects is justified. Therefore the re- 
sults shown in Sec. IIIII all include the proton electric form 
factor with the usual dipole parameterization 



G p E (Q 2 ) 



G D {Q 2 



1 



1 



(12) 



(A = 18.43 fm 2 ). For the neutron electric form factor 
we use the parameterization from [l6[ 



G n E {Q 2 ) = 



Q 



1 + 5.6 



Q 2 w £ 

4m 2 



G|(0 2 ), 



(13) 



with 



-1.911829 fiN and m being the nucleon mass. 



III. RESULTS OF THE LIT CALCULATION 

In this section we present results on i?^, focusing on 
the evolution of dynamical effects as the momentum 
transfer decreases. In Fig. [H we show Rl at constant 
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FIG. 2: (Color online) Longitudinal response function for q = 
200, 150, 100 and 50 MeV/c with the AV18 two-nucleon force 
only (dashed), and with the addition of the UIX (solid) or the 
TM' (dashed-dotted) three-nucleon force. 
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FIG. 3: (Color online) R L (w, q) for 250 < q < 500 MeV/c: calculations with AV18 (dashed) and AV18+UIX (solid). Data 
from [2^] (squares), [2(| (circles), [j| (triangles down), and [l7| (triangles up). 



g = 200 and 100 MeV/c. As already shown in [§] one has 
a large quenching effect due the 3NF, which is strongest 
at lower q. One should notice that such an effect is not 
simply correlated to the under-binding of the AV18 po- 
tential (binding energy Eb= 24.35 MeV in the present 
LIT calculation, with higher EIHH precision E B = 24.27 
MeV 18|). In fact, if this was the case, the results 
with the Malflict-Tjon potential (MT) which gives 
a slight over-binding of 4 He (E B =30.56 MeV), would lay 
even below those obtained with AV18+UIX (E B = 28.40). 
On the contrary the MT curve is situated between the 
curves with and without 3NF. 

We would like to mention that we do not give results 
for the very threshold region, where one has a + res- 



onance [20j. In our present calculation we are not able 
to resolve this narrow resonance, therefore we subtract 
its contribution before inversion as we do for the elastic 
peak 0]. This procedure of course does not affect the 
result above the resonance. 

Given the large 3NF effect at lower q it is interesting 
to see whether there is a dependence of the results on 
the 3NF model itself. To this end we have performed the 
calculation using also the Tucson Melbourne (TM') [2l[ 
three-nucleon force. While the UIX force contains a two- 
pion exchange and a short range phcnomcnological term, 
with two 3NF parameters fitted on the triton binding en- 
ergy and on nuclear matter density (in conjunction with 
the AV18 two-nucleon potential), the TM' force is not ad- 
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justed in this way. It includes two pion exchange terms 
where the coupling constants are taken from pion-nucleon 
scattering data consistently with chiral symmetry. 

Our results with the TM' force are obtained using the 
same model space as for the UIX potential and the ac- 
curacy of the convergence for the LIT is found to be 
at a percentage level, in analogy to that of the UIX as 
described in 0. The cutoff of the TM' force has been 
adjusted on triton binding energy, when used in conjunc- 
tion with the AV18 NN force. With a cutoff A = 4.77771^, 
where is the pion mass, we obtain the following bind- 
ing energies 8.47 MeV ( 3 H) and 28.46 MeV ( 4 He). We 
would like to emphasize that the 4 He binding energy is 
practically the same as for the AV18+UIX case (as al- 
ready found in [22j). 

In addition to what is shown in our recent Letter Q, 
here we investigate also other low-q values. Figure [2] 
shows that the increase of 3NF effects with decreasing q is 
confirmed. Moreover it becomes evident that also the dif- 
ference between the results obtained with two 3NF mod- 
els increases with decreasing q. One actually finds that 
the shift of the peak to higher energies in the case of UIX 
generates for Rl a difference up to about 10% on the left 
hand sides of the peaks. This is a very interesting result. 
It represents the first case of an electromagnetic observ- 
able considerably dependent on the choice of the 3NF. 
In the light of these results it would be very interesting 
to repeat the calculation with EFT two-and three-body 
potentials [HI, [24[ ■ At the same time it would be highly 
desirable to have precise measurements of Rl at low q. 
This could serve either to fix the low-energy constants 
(LEC) of the effective field theory 3NF or to possibly 
discriminate between different nuclear force models. 

In Fig. [31 an overview of the results obtained for larger 
q is given, showing also the comparison with existing 
experimental data. One sees that the 3NF results are 
closer to the data, this is particularly evident at q = 300 
MeV/c. However, the 3NF effect is generally not as large 
as for the lower momentum transfers shown in Figs. [5] 
and [31 In some cases the quenching of the strength 
due to the 3NF is comparable to the size of the er- 
ror bars, particularly for the data from Ref. 25]. The 
largest discrepancies with data are found at q — 500 
MeV/c. While the height of the peak is well reproduced 
by the result with 3NF, the width of the experimental 
peak seems to be somewhat narrower than the theoreti- 
cal one. On the other hand one has to be aware that rel- 
ativistic effects are not completely negligible at q = 500 
MeV/c. They probably play a similar role as found in the 
electro-disintegration of the three-nucleon systems (see 
e.g. Ref. [13). In the case of q = 250 MeV/c the ex- 
perimental results are not sufficiently precise to draw a 
conclusion. 



IV. MULTIPOLE ANALYSIS 

It is interesting to analyze the results of Rl in terms of 
its multipole contributions. Using Eq. ([5]) on the right- 
hand-side of Eq. ([8]), one can decompose £l(v, q) into a 
sum of multipole contributions C J ' x (<r, q) . We have cal- 
culated each of them separately, solving the correspond- 
ing equations ([5]) . After the inversion of the transform we 
have obtained the various multipole responses. By mul- 
tiplying them with the isoscalar/isovector nucleon form 
factors we generate the multipole contributions to the 
longitudinal response function R J £ X (uj,q). In Fig. [?]we 
show how the isoscalar and isovector parts of Rl are 
built up from their multipole contributions at a higher 
(500 MeV/c) and a lower (100 MeV/c) value of q. As ex- 
pected, the higher the momentum transfer, the larger the 
number of multipoles that have to be considered to reach 
convergence. For q — 500 MeV/c up to seven multipoles 
are considered, while for q = 100 MeV/c only three mul- 
tipoles are required for a converged result. Regarding 
the strength distribution among the multipoles two facts 
are evident: i) at higher q the strength is almost equally 
distributed among the first isovector multipoles, while 
in the isoscalar channel the quadrupole gives the largest 
contribution; ii) at low q, as expected, the response is 
dominated by the isovector dipole contribution. While 
the isoscalar dipole is completely negligible, the isoscalar 
quadrupole contributes a few percent. A careful reader 
may notice that the isoscalar response at q = 100 MeV/c 
does not seem to show a convergence in the multipole de- 
composition. As already mentioned, the isoscalar dipole 
is negligible (explaining why the curve labeled with "0" 
is overlapping with the "+1" one), and a similarly neg- 
ligible strength is found for the multipoles higher than 
the quadrupole. This fact is also seen in Fig. \5\ where 
the total strength of the various multipoles is shown for 
q = 100, 300 and 500 MeV/c. At q =100 MeV/c the 
J = 1 and 3 multipoles of the isoscalar response are tiny. 
The J > 3 multipoles are neglected in Figs. [U [2] and Q] for 
the q <100 MeV/c kinematics. We would like to point 
out that the total strengths presented in Fig. [5] do not 
contain the nucleon form factors. The strengths can be 
obtained by integrating in energy (up to infinity) the in- 
version of C J,x (a,q) (non-energy weighted sum rule) or 
just by taking the norm of the right-hand-side of Eq. (JSJ) 
for each multipole, i.e. the norm of C J,X (q)\^o) . 

In Fig. 3NF effects are illustrated in addition. The 
three shown q values give an idea of the evolution of the 
effect from the short range to the long range regime. At 
the highest q value the three strongest contributions are 
given by the isovector dipole and the isoscalar and isovec- 
tor quadrupole. They are enhanced by the 3NF, while all 
other multipoles are decreased, resulting in a net small 
quenching effect. At q= 300 MeV/c one notices a kind 
of transition situation where only the still dominating 
isovector dipole strength is increased by the 3NF, while 
all other multipoles are quenched. At q = 100 MeV/c 
the strength of all multipoles is decreased by the 3NF 
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FIG. 4: Response functions of the lowest isovector (upper panel) and isoscalar (lower panel) Coulomb multipoles, starting 
with the monopole and consecutively adding higher multipoles up to Jmax = 7 for q — 500 MeV/c (left) and J ma x = 2 for 
q = 100 MeV/c (right) in case of the AV18+UIX potential. 



resulting in an overall sizable quenching effect. 

In Fig. [6l the effect described above shows up 
more clearly in the energy distribution of the domi- 
nant multipole contributions (isovector dipole, isoscalar 
quadrupole) at q — 500 and 100 MeV/c. In particular at 
q = 500 MeV/c the increase of the strength due to the 
3NF in the isovector dipole channel is found mainly in the 
high energy tail, while in the isoscalar quadrupole the in- 
crease is found around the peak. At q = 100 MeV/c the 
situation is different in that the quenching of the strength 
due to the 3NF concentrates in the peak region, for both 
multipoles. The net result of this mechanism is the in- 
crease of the 3NF quenching effect with decreasing q that 
was evident in Fig. O 

Here we would like to comment on the fact that the 
large contribution of the 3NF at low q seems at vari- 
ance with the smaller contribution to the photoabsorp- 
tion cross section [18j, dominated as well by the isovec- 
tor dipole. The reason is twofold. It has to do on the 
one hand with the correct use of the Siegert theorem, 
and on the other hand with the common procedure to 
let theoretical cross sections start from the experimental 



threshold, also when the binding energies do not repro- 
duce the experimental values. The detailed explanation 
is in order here. Because of charge conservation the rela- 
tion (Siegert theorem) between the charge dipole matrix 
element (C 1,v ) considered here and the electric dipole 
matrix element (El) considered in the photon case im- 
plies the factor (E n — Eo) (see also [Hj]). The binding en- 
ergy (—Eo) of 4 He, however, is about 15% lower for AV18 
than it is for AV18+UIX. One has the following conse- 
quences. In the AV18+UIX case of Ref. [18j the result 
of the squared matrix element is multiplied by (E n — Eq) 
which is equal to w 7 , while in the AV18 case a multi- 
plication by (E n — Eq) implies a smaller multiplicative 
factor. Therefore the quenching three-body effect results 
to be smaller. It is only thereafter that the AV18 4 He 
total photoabsorption cross section is shifted to the ex- 
perimental threshold. 

One might question the procedure of shifting the the- 
oretical cross sections to the experimental threshold, but 
without such a shift all 3NF effects would be much am- 
plified in the photoabsorption cross section and even in 
the present response function results. However, in this 
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FIG. 5: (Color online) Isovector (left) and isoscalar (right) multipole strength distribution of Rl at g=500, 300 and 100 MeV/c 
in case of the AV18 and AV18+UIX potentials. 



case one could say that they are in a way "trivial binding 
effects". 



rules map an energy dependent observable into a set of 
n discrete values, integral transforms map the same ob- 
servable into a set of continuous values. 



V. INTEGRAL PROPERTIES OF R L 

There are many examples in different fields of physics 
where one is not able to access a certain observable, but 
only some of its integral properties. Sum rules (n th mo- 
ments of the energy distribution) [2^ | are well known ex- 
amples. They contain a certain amount of often very use- 
ful information about the observable, but a limited one. 
The more sum rules one knows, the larger the amount 
of information at disposal. Integral transforms can also 
be viewed as a special form of sum rules. While sum 



Reconstructing the searched observable from its inte- 
gral properties can be very difficult, since very often only 
a limited number of moments are known or in case of inte- 
gral transforms the result of the mapping does not resem- 
ble at all the observable of interest. This is not the case 
for the LIT. An example is illustrated in Fig. [7^ where 
we compare R L (uJ,q = 300MeV/c) for the AV18+UIX 
potential (as in Fig. [3fe) with the corresponding LIT, 
jCl(ct, q), calculated using a typical value of ctj (20 MeV). 
One notices the similarity between the shape of the re- 
sponse function and of its integral transform. This simi- 
larity is due to the fact that the Lorentz kernel is a repre- 
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FIG. 6: (Color online) Response functions of the isovector dipole and isoscalar quadmpole for q — 500 and q — 100 MeV/c 
with the AV18 (dashed) and AV18+UIX (solid) potentials. 



sentation of the (5-function. It is this property that makes 
the inversion of the integral transform [30| reliable and 
sufficiently accurate for the Lorentz kernel. This situ- 
ation has to be confronted with the Laplace transform. 
The Laplace transform of Rl, called Euclidean response, 
is given by [3l[ 



£(r,q) 



du> exp 



g 

2m 



Z\G'W) 



(14) 

Fig. [7b shows how £{t) exhibits a completely different 
form than Rl(co, q)- It is interesting to observe that even 
in the r-space the Euclidean response obtained using only 
the two-body potential gives a result much different from 
that which includes the 3NF. However, what is not evi- 
dent from £(t) is in which energy region the contribution 
of the 3NF is important. 

Fig. [7)3 also shows a comparison of our results with 
those of Ref. [3l| obtained with the Monte Carlo method. 
The comparison is of interest, even if the potentials used 
in the two cases are slightly different (in [3l| the older 
versions of the Argonne and Urbana AV8 and UVIII 
had been used). At r larger than 0.02, our result with 
AV18+UIX lies within the error band of the Monte Carlo 
numerical noise. At smaller r, and in particular at t = 
Fig. [7p shows a discrepancy between the present £{r,q) 
and that of Ref. [3l[ . This is certainly due to the differ- 
ent potentials used. The £ (0) value corresponds to the 
zero-th moment of Ru(u),q). This is a classical integral 



property of Rr,(w,q) that has been much discussed in the 
literature under the name of Coulomb sum rule (CSR) 
(for a review see [2!| HI])- The sum rule consists in con- 
necting the integral of the inelastic longitudinal response 
to the number of protons and to the Fourier transform 
of the proton-proton correlation function p pp (s), i.e the 
probability to find two protons at a distance s. In fact 
for the charge density operator of Eq.([JJ (and neglecting 
the neutron charge form factor) one has 



CSR(g) 



R L {u,q) 



= Z+Z(Z-l)f pp (q)-Z 2 \F(q) 



where uj e i represents the elastic peak energy, f PP (q) is the 
Fourier transform of p pp (s) and F(q) is the nuclear elas- 
tic form factor. (Another interesting sum rule concern- 
ing the second moment of p PP {s) has been considered in 
Ref. [13, where one finds (s 2 ) = 5.67 fm 2 for AV18+UIX) 
The main interest in CSR(g) has been considered its very 
simple model independent large q limit, i.e. the number 
of protons, information about the proton-proton corre- 
lation function. 

In Fig. [51 fvv(q) is shown in comparison with the re- 
sults of Ref. [f| obtained with the same potential, in- 
cluding in addition higher relativistic order corrections 
as well as exchange operators, neglected here. We can 
make the following observations: i) the perfect agree- 
ment of the two results when the operator in Eq. ([JJ is 
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FIG. 7: (Color online) Longitudinal response function Rl 
(solid) at q — 300 MeV/c with its Lorentz Integral Trans- 
form (dashed) for the AV18+UIX potential (a). The LIT has 
been multiplied by ^f- to normalize the integral kernel with 
respect to Eq. (0. Euclidean longitudinal response for the 
same momentum transfer (b): comparison of the GFMC cal- 
culation [3l[ for the AV8+UVIII potential (band between thin 
lines) and the result of this work with the AV18 (dashed) and 
AV18+UIX (solid). 
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FIG. 8: (Color online) (a): f pp (q) in the CSR (Eq. ([15])) 
with AV18 (dashed) and AV18+UIX (solid) potentials. Dots: 
results from [f| with the one-body density operator of Eq.([3jl, 
triangles: results from 0] with two-body density operator, 
(b): percent difference between the curves in (a) 



used shows the great accuracy of the two calculations; 
ii) the contributions of exchange currents, become neg- 
ligible below q = 300 MeV. This means that at low q 
the physical interpretation of f pp (q) as Fourier transform 
of Ppp(s) is safe. Therefore, in principle, the comparison 
theory-experiment would allow to study microscopically 
the largely unknown long range correlations, iii) as one 
can see in Fig.[Sj3 the effect of the 3NF on f pp (q) is up to 
15% in the "safe" region below q=300 MeV/c This gives 
an idea of the required experimental accuracy. 

Unfortunately, obtaining the "experimental" CSR(g) 
(as well as £(Q)) is a non trivial task, due to the ne- 
cessity of extrapolating data up to infinite energies, even 
crossing the photon point, where (e, e') measurements do 
not have access. Different extrapolating functions have 
been proposed. They have been used also recently in 
Ref. [13]. Our results can help to determine these tail 
contributions. They can in fact be obtained subtracting 
from CSR(g) the experimental sum of the data up to the 
last measured point at oj max . From the Saclay data [26| 



at q =300 and 350 MeV/c we have estimated this high en- 
ergy contribution to be about 7% of the CSR. The effect 
becomes twice as large for the higher g-values. However, 
while this procedure would be safe enough at low q, at 
large q this estimate can be inaccurate because of the 
neglect of relativistic effects, two-body operators and the 
role of neutron form factor. In general it would be de- 
sirable that the tail contribution does not overcome the 
3NF effect. Therefore accurate data should be taken as 
far in energy as possible. Of course they cannot overcome 
the photon point, therefore it is interesting to calculate 
the contribution of the tail beyond it. In Table QJ one 
sees that for q- values up to 200 MeV/c the contribution 
of the time-like region remains very low reaching at most 
1.5%. 

The previous discussion gives an idea of the experimen- 
tal accuracy required to access the information about the 
proton-proton correlation function. 
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TABLE I: Theoretical CSR for low q with the 
AV18+UIX potential in comparison with I* h = 
J 1 ^ 9 dLoR e L xp ( qi uj)/\G p E (Q 2 )\ 2 . The percentage contribution 



of the time-like response is also shown. 



q [MeV/c] 


CSR(q) 


rth 
1 1 


% time-like 


50 


3.88 


3.86 


0.5 


100 


3.57 


3.55 


0.6 


150 


3.17 


3.14 


0.9 


200 


2.79 


2.74 


1.5 



VI. CONCLUSIONS 

In this paper we have analyzed 3NF effects on the elec- 
tron scattering longitudinal response function at several 
kinematics. The most interesting results regard momen- 
tum transfers between 50 and 250 MeV/c. Large effects 
of 3NFs are found for two different three-body potentials 
(Urbana IX and Tucson-Melbourne'). We also observe 
that the 3NF effects differ by non negligible amounts for 
the two three-body force models. Since the difference 
between the two results increases with decreasing mo- 
mentum transfer one can ascribe it to rather different 
long range correlations generated by the the two forces. 
This fact results from the difference in the continuum ex- 
citation spectra of the two potentials. In light of this ob- 
servation one can envisage the possibility to discriminate 
also between phenomenological and effective field theory 
potentials, if precise experimental data were available at 



these kinematics. 

Three-body force effects have been analyzed separately 
in the various multipoles contributing to the response. 
Below q = 300 MeV/c a cooperative quenching effect in 
all multipoles has been found. 

Integral properties of the longitudinal response func- 
tion have also been addressed. In particular the possibil- 
ity to extract information about the long range behavior 
of the proton-proton correlation function has been dis- 
cussed. In relation to this it has been underlined how 
the present results can be used in determining the en- 
ergy tail contributions to the Coulomb sum rule. 

In general it has been emphasized that different from 
the search for the short range correlations, the study 
of the long range ones is not affected by complications 
due to relativistic and two-body contributions. There- 
fore a Rosenbluth separation of inclusive electron scatter- 
ing cross section of 4 He at momentum transfer q < 200 
MeV/c would be of high value in view of a more accurate 
determination of the three-body force and in general of 
the long range dynamics of this system. 
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